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1. Introduction. 

In 2003 van Straten and van Enckevort during a computer search found the third order 
differential operator (77) 

L^e^ - x{2e + l)(116l^ + 1161 + 5) + 12hx'^{e + if 

Here 9 — x4^. Then Ly — where 



for some unknown coefficients A„, n = 0,1, 2,. ..It took me five years before I found an 
explicit expression for A„, namely Aq = I and 

/2n\ /3n\ ^ n - 2fc fn\ ^'^ fnXfbn- 5fc\ '^ (5fc)! (5n-5fc)! 

\n) \n) ' ^n-hk\k) \hk)\ n J kl^ {n - k)\^ 

for n > 0. How is it possible to find such a complicated formula? Let us first point out that 
I immediately got an e-mail from Zudilin where he (trivially) simplified the formula to 

k=0 \ / \ 



and later (not so trivially) to 

1/5] / X 3 



fc=0 ^ ^ 



Mn - 5fc - 1\ /An - 5k 
\ 3n J \ 3n 



Here it is evident that An is an integer. But it was found in the complicated form above and 
here we shall tell the long story how I came to let Maple's "Zeilberger" find the recursion 
formula for this monster. 

2. Hypergeometric equations. 

The first Calabi-Yau differential equations connected to Calabi-Yau manifolds needed in 
string theory appeared in the 1980-ies in the physics literature. They were of the form 

- xP{e) 

where P{d) is a polynomial of degree 4. There were 13 of them and I found the 14-th 
equation 

.^-12«.(.-|-1)(.-hA)(.+ ^)(.+ H) 

while writing the popular paper [1]. It was also found independently by C.Doran and 
J.Morgan [14]. It is interesting that the same 14 4-tuples of fractions ( like (1/12, 5/12, 7/12, 11/12)) 
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occur in J.Guillera's Ramanujan-like formulas for I/tt^ (see [14]). The hypergeometric equa- 
tions arc 7^ 1-14 in [3] called the "Big Table" from now on. 
3. Formal definitions. 

After I had finished writing [1] I found the papers [10], [11] by Batyrev, van Straten et al. 
There were another 14 equations 15-28 in [3] ) coming from geometry and they were not 
hypergeometric. What is the common property of these 28 equations? 

Definition: A Calabi-Yau differential equation is a 4-th order differential equation with 
rational coefficients 

+ a3{x)y"' + a2{x)y" + ai{x)y' + ao{x)y = 

satisfying the following conditions. 

1. It is MUM (Maximal Unipotent Monodromy), i.e. the indicial equation at x = has 
zero as a root of order 4. It means that there is a Frobenius solution of the following form 

yo = 1+Aix + A2x'^ + ... 

2/1 = yo log(a;) + Bix + B2x'^ + .. 

2/2 = ^yolog^(a;) + {Bix + B2x'^ + ...) log(a;) + Cix + Caa;^ + ... 

2/3 = ^2/0 log^ {x) + ]-{Bix + B2X^ + ...) log^ {x) + {Cix + C2X^ + ...)\og{x) + Dix + D2X^ + ... 

D 2 

It is very useful that Maple's "formaLsol" produces the four solutions in exactly this form 
(though labelled 1 - 4 ) 

2. The coefficents of the equation satisfy the identity 

11 3 1 

ai = -0203 - -al + a'2- ^0303 - -03 

3. Let t = 2/1/2/0. Then 

q = exp(t) = X + C2X^ + ... 

can be solved 

X = x{q) = q - C2(^ -I- .... 
which is called the "mirror map". We also construct the "Yukawa coupling" defined by 

This can be expanded in a Lambert series 

where the are called "instanton numbers". For small d the rid are conjectured to 
count rational curves of degree d on the corresponding Calabi-Yau manifold. Then the third 
condition is 

(a) j/o has integer coefficients 

(b) q has integer coefRcients 

(c) There is a fixed integer A^o such that all Nqu^ are integers 
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Actually it looks as if conditions 1 and 3a,b would imply conditions 2 and 3c. At least 
we have not found any counter example during six years of search. 



3. Fullbacks of 5-th order equations. 

The condition 2 is equivalent to 



2/0 






yi 


2/2 




y's 




y'l 


y'2 



This means that the six wronskians formed by the four solutions to our Calabi-Yau equation 
reduce to five. Hence they satisfy a 5-th order differential equation 

w^^^ + b4w'^*^ + baw'" + b2w" + biw' + bow = 

The condition 2 for the 4-th order equation leads to a corresponding condition for the 5-th 
order equation 

25 

b2 = %b4 - ^bl + Ib's - §6464 - b'l 



5 25 * 2 ^ 5 

The special equation (coming from number theory, see Zudilin [22] ) 

9^ - Sx{20 + 1)(36'^ + 30 + 1){150^ + 156 + 4) - 3x^(0 + if {30 + 2){30 + 4) 

satisfies this equation and it is possible to find its 4-th order "puUback" (# 32 in the Big 
Table). It has degree 8 since 

0^ + ... + 3^x^{30 + lf{30 + 2f 

with large coefficients. We say that two Calabi-Yau equations arc equivalent if they have the 
same instanton numbers. This is the case under the transformation 



where 



Y{x) = f{x)y{g{x)) 
g{x) = X + a^x^ + ... 



Conversely it is conjectured that equivalence implies the existence of algebraic fix) and 
^(a:) as above (see [4]).Yifan Yang ([21]) suggested a different but equivalent puUback which 
usually cuts the degree in half. It also has a symmetry reducing the number of coefficient 
another 50% (see [6], [7]).. It depends on the following. Let 



Wq = X 



2/0 2/1 

2/0 y'l 



2/0 2/2 



y'o y'2 



u>i = X 

Then we have the following identity, " The double wronskian is almost the square" 



Wo Wi 



= x^yl exp 
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(for a proof see [6]). Thus yo is, up to a factor, equal to the square root of the wronskian 
of wq and wi . We give the Yifan- Yang pullback of # 32 

- X {540(0 + 1)4 + 486(0 + if + || 

+x' |72846(0 + 1)4 + ^{6 + 1)' + J} 
1 14580(6* + ^)* + 12717(6» + + 324| 

+ ^x4(60 + 11)2(60 + 13)2 
16 

The most striking result from the Yifan- Ying pullback are 14 degree 2 equations coming from 
hypergeometric 5-th order equations. We show the following example. Consider 

0^- 4. 12*^.(0 + l)(0+l)(0 + A)(,+ ^)(,+ ii) 

which has the Y-Y pullback 

04 - 144a;(16588804 + 3317760^ + 38649602 + 2206080 + 47711) 

+2223i"a;2(40 + 3)(40 + 4)(60 + 5)(60 + 7) 

an equation not found in the computer search mentioned in the introduction because the 
numbers are too big. 

But unfortunately 5-th order differential equations satisfying 25 are rare (except the 
ones constructed from known 4-th order C-Y equations). In [7] they are listed. In the new 
edition of the Big Table there are two new ones, # 355, 356. 



4. Using Maple for finding and factoring differential equations. 

Using Maple's " Zeilberger" is the best way to find Calabi-Yau differential equations com- 
ing from simple sums of products of binomial coefficients. 
Example: #15. Consider the sum 




Then Maple finds the recursion formula for A„. Here NAn = An+i 
with(SumTools [Hypergeometric] ) ; 
r:=Zeilberger((3n) ! /n! '3*binomial(n,k) ^3,n,k,N) [1] ; 

r := (n + 2)'^N^ - 3 {3(n + 1) + 1)} {3(n + 1) + 2} {7(n + 1)^ + 7{n +1) + 2}N 
-72(3n + l)(3n + 2)(3n + 4)(3n + 5); 

Maple finds the recursion of lowest order (degree in N ), say 

Coin)NP + ci{n)NP-'' + ... + Cp{n) 
Converting to a differential operator we find 

Co{e -p)+xci{0 -p+1) + ... + xPcp{0) 
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Hence the C-Y condition 1 (MUM) is equivalent to 

co(n) = {n + p)^ 

But this is not always the case. 

Example: #22. Consider the sum 




Then we compute 

r:=Zeilberger(binomial(n,k)'5,n,k,N)[l]; 
r co{n)N^ + ci{n)N'^ + C2{n)N + C3(n); 
where 

Coin) = (55n2 + 143n + 94) (n + 3)^ 

This will give a differential equation L of degree 6 which can be factored in Maple. Observe 

that first we have to convert the difi^crcntial equation so that 6 is replaced by This can 
be done in Maple using Stirling numbers (see p. 19 in [6]). 
with(DEtools); 

DFactor(L,[Dx,x], 'one step '); 
L2:=%[2]; 

Here L2 is the right factor of L. 

This strategy usually docs not work if co(n)/(n +p)'* contains only linear factors of the 
form n + a,2n + a,3n + a, 4n + a. There are two exceptions: 
#354 co(n) = (n+l)(7n + 12)(n + 3)^ 
#361 co(n) (2n + 5)(8n + 77)(n + 4)'* 

We have also found three cases where we have an irreducible factor of degree four. 
#251 Coin) = (441n^ + 3780n3 + 11634n2 + 15207n+7202)(n + 4)4 
#299 Coin) = (4704n4 + 4057271^ + 11773271^ + 133251n + 49684)(n + 4)^ 
#367 co(n) = (ISSSn-* - 13592n3 + 35556n2 - 39931n + 16322)(n + 4)" 
There is also a case with an reducible factor of degree 3 
#235 Coin) = (5n + 9)(6n2 + 22n + 19) (n + 4)* 



5. Multiple sums of binomial coefficients. 

There is a " MultiZeilberger" but it is too slow to be of any practical use. So we use 
"brute force" instead. Let us take an 
Example. #349. Let 

..=i:(-i,-3»-(3;)<f(:)(„-)(:) 

Assume that the annihilating differential equation is 

Ly = 502/^"^ + qiy'" + <l2y" + qiy' + 902/ = 

where 

go = a^x'^ + a^x^ + ... + aiQX^^ 



91=63x^ + 640; +... + 615X 
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92 = C2X^ + csa;^ + ... + cmx^* 

13 



qz = dix + d2X +... + dizx 

94 = eo + e\x + ... + 6120;^^ 

with 65 unknown coefBcients 04, ..ei2. Then we compute 70 coefBcients of A„ (takes about 
one second) and substituting 

69 

n=0 

in Ly we get 70 hnear equations which is solved in a few seconds on a laptop. #349 is the 
only known case where we really need coefficients of degree 12. 



6. Hadamard products. 
If 

00 00 
u = V = c„a;" 

n=0 ra=0 

are two D-finite (i.e. satisfying differential equations with polynomial coefficients) power 
series, then the Hadamard product 

00 00 
y = u*v = 'Y^ anX^ = ^ 6„c„a;" 

n=0 n=0 

is also D-finite. It was suggested by Duco van Straten that if u and v satisfied "nice" 
second order equations then u*v could satisfy a Calabi-Yau equation. This was indeed the 
case for the following type 

- x(o6i^ + ae + b) + cx^{e + 1)^ 

with integer a, 6, c, a class of differential equations studied by Don Zagier ([21]). There are 
10 such equations listed as (a),(b),..,(j) in [2] . There are also 10 third order equations 

- x{2e + VfiaO"^ + ae + b) + cx^{e + if 

suitable for Hadamard products with e.g. (^^") giving 5-th order equations. There are 

relations between the coefficients a, 6, c and a,b,c described in [4] There are many 
equivalences between the various Hadamard products. All this treated in [4]. 



7. The mirror at infinity. 

There are many Calabi-Yau equations ending with a term cxP{9 + 1)^. For all these 
equations there is a "mirror at infinity" obtained by the transformation 9 — > —6 — 1 and 
x — > ax~^ for a suitable constant a. E. Rodland did this for # 27 in [19]. 

Example. # 193. We have 
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and the equation 

- 7a;(11356'^ + 22046*^ + 16836*^ + bSlO + 77) 

+a;2 (287236*^ + 4070861^ + li2me^ - 133761 - 896) 
-a;^(321266'^ + 3851461^ + 2651161^ + 1073161 + 1806) 
+7 • lla;''(1306'^ + 2540^ + 19261^ + 656' + 8) + ll2a;^(6i + 1)" 
We make the substitutions 9 — > — ^ — 1 and x — > 2>~'^x~^. It follows 

11^6^ - 7 • lla;(1306'^ + 2666'^ + 21061^ + 776' + 11) 

-a;^(321266'* + 899906'^ + 1037256'^ + 552536' + 11198) 
+a;^ (287236'* + 741846'^ + 634746'^ + 206256' + 1716) 
-7x^(11356'* + 23366'^ + 18816'^ + 7136' + 110) + 7V(6' + 1)* 
which is #198. In this case we also know a formula for the coefficients 




which is rather unusual (11 out of 36 cases). Also equations ending with cxP{29 +1)** can 
be treated in a similar way by the substitutions 6 — > —6 — 1/2 and x — > ax~^. 



8. Harmonic Sums. 

At the end of the paper [18] by P.Paule and C.Schneider there is a remark that using 
"Zeilberger" on 



M„"=5:(n-2fc)Q' 



which is identically zero by symmetry, one obtains the same recursion formula as they ob- 
tained for 

^" = E Q {1 + ^i-^Hk + 7Hn-k)} 

where 

n ^ 

if n > 1 and Hn = if n < 0. The recursion gives the differential equation # 27 in the 
Big Table. 

Lemma 1. We have 

An! = n!(if„-7) 

where 7 is Euler's constant. 
Using this we find 



2dk 



7 



Indeed 

-lit"- ' = s t K:)' = E + . 7.„_,, 

fc=0 ^ ^ fe=0 ^ ^ fe=0 ^ ^ 

In this way 28 equations of type 

"An" =^{n-2k)C{n,k) 

k 

with 

C(n, n — fc) = C(n, k) 
were found, the last being #360 with 

(n + 3fc\ /4n - 3fc\ (3fc)! (3n - 3/c)! 



\ ^ / \ ) k\^ {n — k)\ 
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{1 + k{-AHk + AHr,-k + 3Hn+3k - 3if4„-3fe)} 

But sometimes it is not enough to take the derivative of "An" , we also have to sum over 
negative A;. (this was pointed out to me by Christian Krattenthaler). For this we need 
Lemma 2. Let n be a positive integer. Then 

T{-n + x) = ^-^x-' +0{1) 
ni 

when X 0. 

Proof. We have 



sin(7rt) 

If t = n + 1 — X we obtain 

1 



T{-n + x) = — 



sin(7r(n + 1 — a;)) T{n + 1 — a;) 

= -^^tin^ I = JzlTL + Oil) 

sin(7rx) T{n + 1 - x) xr{n + 1) ^ ' 

To illustrate this we consider 
Example #264. Let 

2rA ^ X fn\ f2k\ f2n - 2k\ (2n + 2k\ ^ [An - 2k\ ^ f2n\ / 2n ^ 



k 

Using Lemma 2 we compute 



n 



k — ej k 



-2k-2e\ 1 f2k\ ^ , _,_2 



2n ^_(-l)V2»+fc^-■^^o(.') 



—k — sj k \ 2n 
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Collecting this we find the derivative of 'M„" at —k and the correct formula 



16" 



En /nW2fcW2n-2feW2n+2fc\2 /4n-2fc\2 /2n\ -1 / 2n \-l 
k=0 \k)\ k ) \ n-k ) \ n+k ) \2n-k) \k) \n-k) 
X {1 + k{ — 2Hk + 2Hn-k — 3Hn+k + ^H2n-k + 2i?2fc — 2il2n-2fc + 4i?2n+2/c — 4if4„_2/c)} 



_l_ Tt+2fc f2n+k\ /2n+2k\ /2n-2k\'^ /in+2k\^ /2k\~ ^ /n+k\~^ (' \ 

~'~ 2^ k \ 2n )\ n+k ) \ n-k ) \ 2n+k ) \k ) In/ \n+k) 



fe=l 



We have a family of equations with 
with 



fe 



7-2p 



2n - 2A:N ^ 
n — k 



p 


s 


# 


#°° 








27 


243 


1 





212 


117 


2 





246 


247 


3 





-6* 




4 





~6* 




5 


2 


247 


246 


6 


8 


117 


212 


7 


14 


243 


27 



For #117 the technique used above gives 

-5/2fc\6/2ri-2fe\'^ 



A = 256- n=o©"°a'^) e:::f) {l + fc(-TJ?. + 7Jy„_. + 12J7,, 



12iJ. 



2n-2fe 



)} 



which is nonsense since wc cannot sum the infinite sum (maybe PARI can do it?). Indeed 

let n = 1 in the second sum. Then 



k=l 



2k 
1 + k 



57 1 



log(2))(a;-l)+0((a;-l)2 



where the constant term exactly cancels the value of the first sum. The same occurs when 
n = 2 and n = 3. So there is still no formula known for An for # 117. There are 12 
other cases (like #243) where taking the derivative with respect to k does not work, but 
for which there exist other formulas for An- 

9, Empty sums. 

Example # 133. Consider the sum 



"4 " — 



E(n-2fc) 



n \ (3A;)! (3n-3fc)! 
3n-3A;y k\^ (n-fc)!^ 



" Zeilberger" gives the recursion 



(n + 2)'^N^ - 12(2n + 3)^(371^ + 9n + 7)N + 432(2n + lf{2n + 3)^ 
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which corresponds to the differential equation 

6*^ - 12x{2e + 1)2(36'^ + 36* + 1) + 432x^(26' + 1)2(26* + 3)^ 

which we recognize as the Hadamard product A* f (see [2]). 

Looking at "An" we observe that in order to (^) to be nonzero we need k < n/3. 
Similarly {^n-3k) nonzero only for k > 2n/3. So the sum is not only zero, it is also 
empty. Consider the case k < n/3. Then 



3n — 3fc — 3e 



7 I 



r(3n - 3fc - 3£ + l)r(-(2n - 3fc - 1) - Se) 

2n-3fc-l 1 i\k /n nj\ -1 



(_l)2n-3fc-l ^ ^ 3(-l)^ /3n - 3k 

{3n - 3k)r {2n - 3k - 1)\ {-3ey ~2n-3k\ n 



which gives 



Ar. 



2 [n/3] 



E(-i)' 



which simplifies to 



fe=o 



An — 3 



n — 2k fn 
2n-3k\k 



3n - 3k\ ^ (3fc)! (3n - 3fc)! 
n ) 'W {n-k)\^ 



' ^ ' .f,n — 2kl'n\^f2n — 3k 
' 2n - 3k U 

fe=0 ^ 



for n>0. 

This is just a special case of 



"^""=^"B-2^)(:)V/=)(3;)( 



n \ {3ky. {3n - 3A:)! 
y3n - 3k) {n-k)\^ 



with the following table 



# 


a 


C'n 


E{n,k) 


133 


-1 




1 


279 


1 




1 


334 


1 




(2n\-^ / 2n \ 
\k ) \n-k) 


281 


2 




(2n\-'^ / 2n \ 
\k ) \n-k) 


363 


-3 




/2n\ / 2n \ 
\ k ) \n-k) 


352 


-1 




f2k\ r2n-2k\ 
\k)\ n-k ) 


253 







(2k\ (2n-2k\ (2n\-^ 
\k)\ n-k J \2k) 


353 


-1 




(2k-\ (2n-2k\ (2n\-^ 
\k) \ n-k )\2k) 


350 


-2 


c:f 


(2k\ (2n-2k\ (n+k\~^ (2n-k\ ~^ 
\k) \ n-k )\ n I \ n I 



The computations before show that 



[n/31 



A„ = 3C„^(-1)' 



fc=0 



n — 2k fn 
2n-3k\k 



3+a 



E{n, k) 



2n - 3k 
n 
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for n > 1. 

Replacing 3 by 4 we get 



" A " 
with the table 



■Ak) k\^ {n-kyA 



# 


a 




E{n,k) 


300 


-2 


(5n)! 




(2r))!r)!-' 


3G 


-2 


o 


1 


364 


-1 




1 


357 







1 


205 





1 


\2kJ 


365 


-3 


1 


r2k\ rin-'2,k\ 
\kl\ n-k 1 



Using Lemma 2 one easily shows 

f2n 



An = 4. 

Replacing 4 by 5 we get 

'M„" =C„^(n-2A:) 

k 

with the table 



[n/4] 
fe=0 



n-2k / 
^ 3n - 4k U 

— n ^ 



4+0 



'3n - 4A; 
2n 



E{n,k) 



n \ {5k)\ {5n - 5ky. 
5k) [^5n-5k) k\^ {n-k)\^ 



# 


a 




E{n,k) 


354 


-3 


1 


1 


B *r] 


-2 


1 


1 



We get 



^-C:)C„>.'|'(-i)'i;^(:)"°.(n,*,( 



4n — 5fc 
3n 



for n > 0. For the case B * rj we get since B* is multiplication by = (^j^) (^^) the 

formula for the coefficient for the third order equation r] 



^„ = 5^(-l)^--- 



fe=0 



4n — 5k\k 



4n — 5fc 
3n 



the equation mentioned in the Introduction. 
Finally there is the case #347 with 



A„ = 6 



n — 2k In 



E 

fc=0 



5n — 6fc V fc 



5n — 6A; 
3n 



One of the most intricate cases is #305 with 



n + 2A;^ (2,n - 2k\ {2n + Ak\ fQn - 4fc 
k 



3n 

n-A; )\n + 2k )\Zn-2k) \n + k 
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There one has to have different sums for —k depending on if k <n/2 or fc > n/2. See the 
final result in its full glory in the Big Table. 



10. Reflexive polytopes. 

M.Kreuzer and H.Skarkc have classified reflexive polyhedra in four dimensions . They 
found 473 800 652 of them. For each of them is associated a Laurent polynomial S in 
four variables Then a solution to a Calabi-Yau differential equation is constructed with 
coefficients A„ =constant term(5"). 

We will show the idea with an example in dimension two where there are only 16 reflexive 
polytopes 

Example. Consider the polygon with four vertices 



with associated Laurent polynom 



2 0-1-1 
-110 -1 



a 11 

S= — +y+- 



xy 



Then 



— (.")=x;(3;)(t)(„!g 



Using "Zeilberger" we get a differential equation of order six which factors into a huge left 
factor of order four and a right factor 

561^ + xe{lie - 1) + 6x^e^ + x^e{13e - 9) - a;"* (29861^ + 16366' + 960) 
-4a;^(7266i^ + 351061 + 2391) - 8a;*^(9176»^ + 475261 + 3519) 
-32a;^(2566»^ + 137261 + 1077) - 32x^(32761^ + 8106' + 528) 
-576x^(5461^ + 7861 + 29) - 128a;^°(6i + 1)(3556' + 417) 
-23808x"(6' + l)((9 + 2) 



with solution 



t/o = 1 + 12x* + 60x^ + 420x^ + ... 



This example shows that a very simple polj^ope can give a differential equation of very high 
degree. 

In dimension four we consider 

Example #325. Consider the polytope with 13 vertices giving the Laurent polynomial 

„ 1 y X z X yz x t ,^ ^ x , 1 1 , 

S=- + ^ + - + - + - + ^ + — + -{l + y + yz) + -il + - + —) 
xxyxzxyzx t y yz 

Then there are only even terms so we compute 

An = c.t.(52") 
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A direct approach by expanding the powers of S and then take the constant term costs a 
lot of computer time. Here we will eliminate t and thus reducing the computer time by a 
factor 500-1000. Let 

1 y X z X yz x 
«:=- + - + - + - + - + — + — 

X X y X z X yz 

p = -{i + y + yz) 

X 

a; 1 1 , 
g= - 1 + - + — 
t y yz 

V = pq = {1 + y + yz){l + - + —) 

y yz 



Then we have 



i+j+k=2n ■' 



To get rid of t we need j = k.lt results 

(2n)! ,^ (2n)! 

-^u (p9)^)=c.t.(_^- 

i+2j=2n n=0 ■ 

To find the equation for 

e 

we need 30 coefficients which are computed in about five minutes on a laptop. Arne Meurman 
has computed an explicit formula for A„ summing over eight indices with complicated 
summation limits, making it not very useful. The equation is 

19^61^ - 19a;(43336i'' + 621261^ + 477861^ + 167261 + 228) 

+a;2 (43074956'^ + 760048461^ + 621640661^ + 2802424^ + 530556) 
-a;^ (9373936961* + 213316800(9^ + 23603719661^ + 12574861261+25260804) 
+a;*(2408138006l* + 7785292006*^ + 1041447759^ + 63180280961+ 138510993) 
-2^ • 409a;^(e + l)(28513246l^ + 100355^^ + 112212416* + 3481470) 
+2^ • 3^ ■ 19^ ■ 409^a;^(6l + 1)(6' + 2)(26l + 1)(26» + 5) 



Let us consider another example,with 18 vertices, vl8. 16766 in the notation of Batyrev- 
Kreuzer [12] with 

S = x{l + t + yt + zt + yzt^) + -{l + \ ^ ^ 



X t yzt"^ 

111 1 1 

+y+-+z+- + -+yt+ — + zt+ — +yzt 
y z t yt zt 

Kreuzer computed 60 coefficients with brute force (time 760 hours). Using the the method 
above eliminating x it took only 71 minutes. Finally Duco van Straten, Pavel Metelitsyn 
and Elmar Schomer using modular arithmetic computed 272 coefficients. They succeded in 
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finding a differential equation of order 6 and degree 25. It is not MUM and does not factor. 
This example leaves some doubts about the reflexive polytopcs. 

So far we know polytopes giving the coefficients of the following equations in the Big 
Table (see 

1-14, 16, 24, 25, 26, 29, 42, 51, 70, ~101, 185, 188, 206, 209, 214, 218, 287, 308, 309, 
324,325, 326, 327 

Only the last four are new. 



11. Other equations. 
11.1. Bessel moments. 

Consider the Bessel moments (see D.Bailey, J.Borwein, D.Broadhurst and M.L.Glasser 

[9]) 

Cm,k = j x^Ko{x)'^dx 



Here Kq{x) is a certain Bessel function that conveniently can be defined by 

oo 

Ko{x)= J e-^"°"'^(*)dt 



This leads to another representation (in Ising theory) 



oo oo 



k\ f f dxi...dx. 



2"* y "J (cosh(a;i) + ... + cosh(a;TO))*^+i 





Let 

15^" 

and 

oo 

Then Yq satisfies the differential equation 

e'^{e - if - 4x6^(2596'^ + 26) + 3600a;^(356»'^ + 706*^ + 636»^ + 2861 + 5) - 3240000x^(6* + 1)* 

The shape of the last term suggests that converting to a; = oo could give a Calabi-Yau 
equation. Indeed 6 — > —9 — 1 and x — > 900x~^ gives the equation 

61^ - a;(356»* + 706»^ + 6361^ + 286» + 5) 

+x^{e + 1)2(25961^ + 5186» + 285) - 225x^(61 + 1)^(61 + 2f, 

e 

which we recognize as # 34 (found by H.Verrill [19]) with solution yo = AnX" with 



n=0 



i+j+k+l+m—n 



n< 



i\j\kW.m\ 
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Similarly 

dn = r5-C6,2n-l 

leads to a differential equation whose mirror at infinity is the 5-th order equation #130 (also 
found by Verrill) with coefficients 



E 

i+j+k-\-l-\-m+s=n 



i\j\kW.m\s\ 



11.2. Differential equations coming from combinatorics. 

There are probably many Calabi-Yau equations whose coefficients count something in 
combinatorics. The last equation # 366 in the Big Table comes from counting random walks 
in Z** (communicated to me by T.Guttmann, [17] ) 

+ xe{39e^ - soe^ - m - 4) 

+2a;2(166'^ - 10706*^ - 105761^ - 6760 - 192) 

-2'^S'^x^{W + 2)(1716'3 + 5666|2 + 6006* + 316) 

-2^3^(3846'^ + 154261^ + 26356^^ + 21736* + 702) 

-2^3^x^{0 + 1)(13936'^ + 55716*^ + 83786* + 4584) 

-2^°3^a;^(6' + 1)(6' + 2)(316'2 + 1056^ + 98) 

-2i23V(6' + l)(6i + 2)2(61 + 3) 

It was found by computing 40 coefficients. An explicit formula for the coefficients is unknown. 
Also # 16 comes from combinatorics, see [16]. 



12. Some remarks. 

Let p be a prime. Expand in base p 

n = no+ nip + n2P^ + ... 
k = ko + kip + k2P^ + ... 

Then it is wellknown that 

Kira Samol and Duco van Straten have found that a similar congruences (Dwork congruences) 
are valid for the coefficients A„ of most Calabi-Yau equations, namely 

An = AnoAniAn^...mod p 

It seems to be valid also for very complicated coefficients, like #264, 274 and also for # 
366 where we have no formula for A„. This property is not preserved under equivalence 
transformations so it is not surprising that it is not valid for some pullbacks of fifth order 
equations. 
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Hundreds of binomial identities resulted from the search of Calabi-Yau differential equa- 
tions. We give only a small sample of the simplest identities 
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